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We report the spin Knight shift {Ks) and the nuclear spin-lattice relaxation rate (1/Ti) in the vortex state as a function of 
magnetic field {H) up to 28 T in the high-Tc superconductor TlSr2CaCu206.8 (Tc=68 K). At low temperatures well below Tc, 
both Ks and l/Ti measured around the middle point between two nearest vortices (saddle point) increase substantially with 
increasing field, which indicate that the quasiparticle states with an ungapped spectrum are extended outside the vortex cores 
in a d-wave superconductor. The density of states (DOS) around the saddle point is found to be kAq^/ H/ Hc2, with k=0.5~0.7 
and No being the normal-state DOS. 



PACS: 74.25.Jb, 74.60.Ec, 74.72.Fq, 76.60.-k 



The vortex state in unconventional superconductors, 
such as the high transition-temperature (high-Tc) super- 
conducting copper oxides which have a c?2.2_j,2-wave gap 
symmetry is of great interest. Among various issues, 
the quasiparticle state associated with the vortices still 
remains far from being well understood. In an s-wave 
superconductor with isotropic gap, the low-lying quasi- 
particle states are localized within the vortex core whose 
radius is ^ = wf/ttAo, with energies E — /iAq/Sf ( 
/i=l/2,3/2....) 1^, where Aq is the energy gap, vp and 
Ep are the velocity and the energy level at the Fermi 
surface, respectively. This prediction was verified exper- 
imentally by the scanning tunneling microscopy (STM) 
1^. In contrast, the gap for d-wave superconductivity 
vanishes in the nodes which leads to a divergent ^ along 
the nodal directions. The corresponding quasiparticle 
state is, therefore, expected to have very different proper- 
ties. Recent theories suggest that the quasiparticle states 
may be extended far outside the cores along the nodal di- 
rections |§,||,^j7|j^. However, this prediction has not yet 
been unambiguously confirmed by experiments. Heat ca- 
pacity measurement at low fields in the vortex state of 
YBa2Cu307 (YBCO) has revealed a ^/H dependence of 
the specific heat coefficient which was taken as evi- 
dence for supporting the theoretical prediction Q|. But 
heat capacity also contains contributions from the vor- 
tex cores which may have a -s/ff-like variation in the 
low field regime as observed in s-wave superconductors 
and pointed out by Sonier et al Our earlier NMR 

measurement on an underdoped cuprate |Tl| and thermal 
transport measurement in the low H/Hc2 regimes are 
consistent with the existence of finite density of quasipar- 
ticle states outside the cores, but no quantitative estima- 
tion of the field-induced states has been available. The 
most serious question about the quasiparticle picture in 
a d-wave superconductor is raised recently by STM ex- 
periments that find no quasiparticle states outside the 



cores ||T^,|T^. Clearly, more spatially-resolved and bulk- 
sensitive studies, and also at higher fields, are needed to 
resolve this issue. 

We report in this Letter measurements of the spin sus- 
ceptibility and the nuclear spin-lattice relaxation rate 
(l/Ti) outside the vortex cores as a function of H up to 
H^28 T (or iJ/iJc2=0.65) in a d-wave high-T^ cuprate. 
We find compelling evidence for the quasiparticles in the 
vortex state being delocalized. Unlike specific heat mea- 
surement that picks up contributions from quasiparticles 
both inside and outside the vortex cores, NMR technique 
can measure separately the quasiparticle states in- and 
out-side the vortex core. In the mixed state, those nu- 
clei located in between two nearest vortex cores have the 
largest distribution probability that gives rise to a singu- 
larity in the field distribution function, while those in the 
cores contribute to a " knee" that is far away from the sin- 
gularity p^ ]. In practice, however, both the singularity 
and the knee are usually smeared out owing to, e.g. lat- 
tice disorders, resulting in a more symmetric NMR line 
shape. Nevertheless, some information about the posi- 
tion selectiveness is retained; the peak of the NMR line 
that results from smeared singularity corresponds to the 
nuclei located outside the vortex cores and in the super- 
conducting state. 

The sample chosen for this study is a slightly over- 
doped cuprate TlSr2CaCu206.8 (T'c=68 K). This com- 
pound has several advantages over the most-extensively- 
studied YBCO family with regards to our purpose. First, 
the Knight shift arising from the spin susceptibility when 
field is applied parallel to c-axis is large , which allows 
an accurate measurement of its change in the supercon- 
ducting state in the presence of magnetic fields. This is 
to be compared to the YBCO family where the corre- 
sponding shift is negligibly small because of the acciden- 
tal cancellation of the hyperfine field, so that measure- 
ment along this direction is not possible Second, 



1 



the spin susceptibility in the normal state above Tc is al- 
most temperature independent, which helps to evaluate 
the density of states (DOS) in the normal state. Third, 
the upper critical field Hc2 ~ 43 T is relatively small , 
so that the field-induced DOS at a given H, if any, should 
be larger. Finally, there is only one crystallographic site 
of Cu, which greatly simplifies NMR spectroscopy and 
enhances data accuracy. 

The polycrystals were c-axis aligned under a field of 16 
T and fixed with a polymer. For all measurements, the 
external field is applied along the c-axis. Figure 1 shows 
a typical NMR line shape at T—A.2 K (gray curve). The 
dotted curve is the corresponding theoretical field distri- 
bution for 75° vortex lattice (VL) obtained as described 
below. The position-dependent field H{r) in the sample 
was calculated by using the London model |17|, 



H{r)=Hj2- 



exp(-GLCV2)exp(-zG, 



1 + 



(1) 



where H is the averaged (externally applied) field, ^ and 
A are the coherence length and penetration depth, re- 
spectively. The summation runs over all reciprocal VLs 

= 2^^SZ(mi +^{l-m- Cosm, where /3 



is the angle between two primitive VL vectors and x and 
y are the unit vectors of the reciprocal VL. We adopt 
P = 75° = 57r/12, which is the value found in neutron 
scattering ^ and STM ^ experiments. The NMR 
spectrum is the spatial average of a Dirac delta func- 
tion, f{H') = {S{H' - H{r)))r. By use of the relation of 
Hc2 = '/'o/27r^^, eq. (1) is left with a single parameter, 
A, that is to be determined by experiment. In the theo- 
retical curve of Fig. 1 (dotted curve) calculated from eq. 
(1), the singularity S corresponds to the middle point 
between two nearest vortices (saddle point), M and m 
correspond to the vortex core and the middle point be- 
tween four vortices, respectively. The solid thick curve is 
a calculation by convoluting eq. (1) with a Lorentzian 
broadening function p{H') = -^rj^^^^Tjjn , with a fixed 
value of (T/7r=170 Oe that is the full width at half maxi- 
mum (FWHM) of the NMR spectrum at r=70 K (above 
Tc). In performing the fitting, both the experimental 
spectrum and the theoretical field distribution curve are 
normalized to unity. From the peak position, we deter- 
mined the Knight shift, Kc^ around the saddle point. We 
also obtained A=950 ± lOOl at r=4.2 K, and hence the 
diamagnetic field Hdia — Hsaddie ~ H that depends on 
the magnitude of externally applied field. Similar value 
of A was obtained from the analysis of ^os^pj lineshape 
which has a narrower width, with FWHM = 50 Oe at 
T=70 K and 140 Oe at T=4.2 K, respectively. 

The value of A=950±100 A is slightly smaller than 
that estimated for the YBCO family (1200-1400 1) 
This is consistent with the current compound being over- 
doped. Also, we have confirmed that the spatially depen- 



dent information is indeed retained, even though the sin- 
gularity and the knee are smeared out. This was done by 
measuring l/Ti in different positions. l/Ti is larger by a 
factor of two in the higher end of the spectrum shown in 
Fig. 1 around the knee, than that measured at the peak. 

Now we turn to the main result, the H dependence 
of the Knight shift and 1/Ti measured at the spectrum 
peak. Figure 2 shows the typical data sets for the T and 
H dependencies of Kc- In the normal state Kc shows very 
weak T-dependence and does not depend on H above 
T=85 K. In particular, Kc is constant (=1.42% ) in the 
temperature range of 85K < T < 150 K. In contrast, at 
low temperatures below the shift strongly depends on 
H. This is clearly shown in Fig. 3 where data at T=A.2 
K are plotted. The solid circles show the data corrected 
by the diamagnetic shift Hdia/H obtained as described 
earlier. It can be seen clearly that Kc — Kdia increases 
with increasing field. This corrected shift is composed of 
the spin part {Kg) and the part due to orbital suscepti- 
bility Korb, respectively, the latter being i/-independent 
in the range of field studied. Note that Kg is related 
to the spin susceptibility and hence the DOS, N{E) as 
Ks(x J N{E){-^§§^)dE, where f{E) is the Fermi dis- 
tribution function. In particular, Kg at low T is propor- 
tional to the DOS at the Fermi level, iV(0). Thus, our 
result indicates that magnetic field induces a finite N{0), 
in the superconducting state. 

This conclusion is corroborated by the TJ-induced in- 
crease of 1 /Ti measured at the spectrum peak. Figure 4 
shows the typical data sets of 1/Ti under various fields 
from H=0 to 28 T. The measurement was done by us- 
ing a single saturation pulse and at the NQR spec- 
trum peak of 22.5 MHz for H=0. As seen in the fig- 
ure, 1/Ti becomes proportional to T below T — 10 K, 
and its magnitude increases with increasing H. This 
is seen more clearly in Fig. 5 where 1/TiT at T=4.2 
K is displayed. In terms of DOS, Ti is expressed as 
l/Ti = a(T) J N{Eff{E){l - f{E))dE, where a{T) is 
the enhancement factor due to antiferromagnetic (AF) 
correlation. For a d-wave superconductor [Q, a{T) was 
shown to have a very weak T variation below Tc and 
q(0) - 0.9a(rc). Therefore, l/TiT oc N(Qf at low T; 
the Ti result also indicates that H induces a finite DOS 
outside the vortex cores. 

Next we evaluate more quantitatively the DOS around 
the saddle point, 5Nh- First, one notices from Fig. 5 
that at high fields of H >5 T, 1/Ti is linear in H which 
suggests that the induced DOS is proportional to V^. 
In fact, in this field regime, a functional of 



5Nh = nN^^BjKk 



(2) 



can fit both the Knight shift and 1/Ti quite well, where 
iVo is the DOS in the normal state. In Fig. 3, the 
solid curve is a fit to Kc — Kdia ~ + Kh\/ H/ Hc2, 
with resultant fitting parameters Kq — 1.00 ±0.03% and 
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Kh = 0.32 ± 0.05%. The solid line in Fig. 4 is = 
0.018 + 3.175-£j [Sec-^K-i]. On the other hand, 1/TiT 
remains finite at H=Q and rises more slowly with increas- 
ing H at low fields below 5 T. We ascribe the zero- field 
residual 1/TiT to that coming from the residual DOS, 
SNimp, due to scattering by impurity or crystal imperfec- 
tion, which is estimated as '^^r""' = \ „ ^^"^J"" — =0.14. 

' JVo V "■^(-'1-' )4.2K 

Then the slower rise of 1/TiT below 5 T can be naturally 
understood as vortex flux preferentially occupying the 
impurity sites where a finite DOS has already been pro- 
duced |21| . In the impurity-dominated regime, a recent 
calculation shows that the field-induced DOS has an 
asymptotic form of 5Nh oc -^-^ln{^j-) where Ha is a con- 
stant (- ^^^K-i). Indeed, the data of <5 T can be 
fitted by such functional as seen in Fig. 5, where the 
resulting dotted curve is ^^r = 0.229 1.650(-£^/7t.^)2 
in unit of Sec~^K~^. The same 5Nh functional fitting 
to the shift of <5 T with Fo=270 T predicts, as 
shown by the dotted curve in Fig. 3, the shift to remain 
at 1.03% when H—Q. Since this value is a sum of the 
impurity-induced spin shift and Korb, namely, 1.03% = 
^^^{1A2%- Korb) + Korb, we obtain Korb=0-96% and 
hence i^s=0.46%, with uncertainties of ±0.03% (abso- 
lute value) Then, k = Kh/Ks is estimated to be 
0.70±0.06%. On the other hand, from the Ti result and 
following the same procedure of estimating SNimp, the 
prefactor k is estimated to be 0.50±0.03%. This differ- 
ence of K by 30% might arise from the choosing of the AF 
enhancement factor a(0) ~ 0.9a(rc) that underestimates 
SNh from the Ti result. The analysis of the Knight shift 
data could also overestimate k; if a triangular model is 
assumed for the VLs, k would be reduced by ^0.1. 

To summarize the above discussions, it is important to 
stress that both the if-enhanced Knight shift and l/Ti 
observed around the saddle point can be consistently ac- 
counted for by a DOS of SNh — nN^yJ H / Hc2, with a 
conservative estimate of/t = 0.5~0.7. We therefore con- 
clude that the quasiparticle states are extended outside 
the vortex cores and that they have a continuous exci- 
tation spectrum without a gap. This feature is in sharp 
contrast with that in an s-wave superconductor. Volovik 
first pointed out that the spatially averaged DOS associ- 
ated with a single vortex in a d-wave superconductor is 
SNh — kNq H/Hc2, where k is on the order of unity [|| . 
Although no calculation of the DOS outside the vortex 
cores is available to allow a direct comparison to our re- 
sults, the substantially smaller k, than unity found in our 
experiment could be due to the effects of vortex lattice, 
its disorder | |2^ , or due to induction of second component 
of the order parameter |^ . 

In conclusion, from l/Ti and the spin Knight shift 
measured around the saddle point in a slightly overdoped 
high- Tc cuprate, we find that the quasiparticle states are 
extended outside the vortex cores with ungapped spec- 



trum in a d-wave superconductor, which is in contrast to 
the localized states within the vortex cores in s-wave su- 
perconductors. The density of such quasiparticle states 
near the saddle point, in the field range of 5 T< H <28 T, 
follows kNq^J H I Hc2 where Nq is the normal-state DOS. 
The prefactor k is estimated to be in the range of 0.5^0.7. 
Although such states have not been found in the STM 
experiments, one may be able to reconcile the STM and 
NMR results by considering the momentum-dependent 
tunneling matrix in the STM experiments which likely 
vanishes in the nodal directions [ p6[ . 
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FIG. 1. ''^Cu NMR line shape at H= 4 T in the vortex 
state (gray curve) where the horizontal axis is measured from 
/=45.14 MHz. The dotted curve is the calculated field dis- 
tribution for 75° vortex lattice. The thick solid curve is a 
convolution of the dotted one with a Lorentzian broadening 
function (see text). 

FIG. 2. Temperature and magnetic field dependencies of 
the ^''Cu Knight shift. The arrow indicates the supercon- 
ducting transition temperature at zero field, Tcq. 



FIG. 3. Magnetic field dependence of the '''^Cu Knight shift 
at T=4.2 K. The open and solid circles are the raw data and 
those after correction for the diamagnetic shift, respectively. 
The error bars arise from the systematic error in estimat- 
ing A. The solid curve is a fit of the data above 5 T to 
Kc{4..2K) - Kd^a = Ko + Kh^JH/Hc2- The broken curve 
is expected when impurity-scattering dominates (see text). 



FIG. 4. Temperature and magnetic field dependencies of 
l/Ti. Note that l/Ti oc T, as indicated by the broken line, 
below 10 K at all fields. 



FIG. 5. H dependence of l/TiT at r=4.2 K taken from 
Fig. 4. The solid line is a linear fit of the data of H >5 T. The 
broken curve is j^r = 0-229 + 1.650(-^Zn2^)^ [Sec'^K"^]. 
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